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ABSTRACT 

We present the abundance analysis for a sample of 7 red giant branch stars in the metal-poor globular cluster NGC 4372 based 
on UVES spectra acquired as part of the Gaia-ESO Survey. This is the first extensive study of this cluster from high resolution 
spectroscopy. We derive abundances of O, Na, Mg, Al, Si, Ca, Sc, Ti, Ee, Cr, Ni, Y, Ba, and La. We find a metallicity of [Ee/H] = 
-2.19 ± 0.03 and find no evidence for a metallicity spread. This metallicity makes NGC 4372 one of the most metal-poor galactic 
globular clusters. We also find an a-enhancement typical of halo globular clusters at this metallicity. Significant spreads are observed 
in the abundances of light elements. In particular we find a Na-0 anti-correlation. Abundances of O are relatively high compared with 
other globular clusters. This could indicate that NGC 4372 was formed in an environment with high O for its metallicity. A Mg-Al 
spread is also present which spans a range of more than 0.5 dex in Al abundances. Na is correlated with Al and Mg abundances at a 
lower significance level. This pattern suggests that the Mg-Al burning cycle is active. This behavior can also be seen in giant stars of 
other massive, metal-poor clusters. A relation between light and heavy ^-process elements has been identified. 
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1. Introduction 

Globular clusters (GCs) provide an exceptional laboratory for 
studying the star formation history of a galaxy because they are 
powerful tracers of the various components (halo, thick disk and 
bulge). Moreover, the ages, abundances and kinematics of these 
objects represent the fossil record of the galaxy formation pro¬ 

* e-mail: isanroman@cefca.es 


cess, since they are amongst the oldest objects in the Universe. 
Starting several decades ago, spectroscopy as well as photome- 
try (e.g. iPiotto et ^l2007t Gratton, Carretta & Bragagli^l2012h 
have shown that GCs are much more complex than previously 
imagined, in particular with the almost ubiquitous discovery of 
multiple populations. None of the scenarios proposed to date 
can fully account for the abundance trends observed in Galac¬ 
tic Globular Clusters (GGCs). Detailed chemical studies of large 


Article number, page 1 of 1 141 







A&A proofs: manuscript no. Sanroman_NGC4372_accepted 


samples of GGCs are required to uncover the evolution of these 
objects. Providing high quality data on more GCs is clearly 
needed to strengthen any conclusion. We present here results of 
the first ever high-resolution spectroscopic analysis of the globu¬ 
lar cluster NG C 4372 based on observations from the Gaia-ES O 
survey (GES, iGilmore et alJl201^[Randich & Gilmorell2012h . 

The GES is a public spectroscopic survey using the high- 
resolution multi-object spectrograph ELAMES on the Very 
Large Telescope. Targeting more than 10^ stars, GES covers all 
the major components of the Milky Way, from halo to star form¬ 
ing regions, providing a homogeneous overview of the distribu¬ 
tion of kinematics and elemental abundances. The survey also 
maps a significant sample of more than 80 open clusters, cover¬ 
ing all the accessible cluster ages and masses. 

Ensuring that GES has maximal legacy impact is of key pri¬ 
ority. Thus the survey has identified a suitable set of objects and 
fields for the calibration of the data in terms of astrophysical pa¬ 
rameters and abundance ratios (Pancino et al. 2012). By includ¬ 
ing objects and fields that are a) well studied in the literature; 

b) in common with other large spectroscopic surveys; and/or 

c) with extremely well measured properties; the results can be 
put in context, compared with other surveys and when available 
combine different data-sets. The GES includes GGCs as cali¬ 
brators, adding some relatively unstudied clusters to the ones 
already present in the ESO archive. To maximize the scientific 
output from the calibrator data-set, a good balance of mandatory, 
known and new/interesting GGCs has been included in the target 
selection. Within this context, this study is focused on the chem¬ 
ical abundances of NGC 4372. With limited spectroscopic data, 
nothing is known about any possible abundance variations of this 
metal-poor halo cluster. We present here an abundance analysis 
of a large number of elements (eg. light elements, a-elements, 
iron peak elements and neutron-capture elements) and analyze 
its stellar population. 

The vast majority of old GGCs stu died in detail to date , 
with the possible exception of Rup 106 (IVillanova et al.llT013h . 
show the chemical signatures of hosting (at le ast) two stel - 
lar po pulations. Several other cases like Pal 12 (ICohenI 120041) . 
Ter 7 (iTa utvaisie ne et al.l 120041 ; ISbordone et al.l l2005h or Ter 8 
(ICarretta et alJ201^ have been proposed as potential single pop¬ 
ulation globular clusters but the small sample investigated does 
not allow for clear conclusions. We also note that all these ob¬ 
jects except for Rup 106 are associated to the Sgr dwarf galaxy. 
There has been much improvement on the study of multiple 
populations within GGCs. The main evidence for this com¬ 
plexity comes from the presence of chemical inhomogeneities 
that in most cases are limited to light elements. Light elements 
like Li, C, N, O, Na, Mg, Al, are known to (anti-)correlate. 
The most outstanding signature is the Na-O anti-correl ation, de¬ 
tected so far in every GC study (ICarretta et al.ll20()9alf3l . again 
with the likely exception of Rup 106. These anomalies have 
been observed also in old, massive extra-galactic GCs in Eornax 
(Letartee^_^ 20061) and in the Large Magellanic Cloud (EMC) 
( Mucciarelli et al.ll2009l). but not in intermediate-age EMC clus¬ 
ters (iMucciarelli et alJ|2014l) . A similar feature, often observed 
in s ome GCs but not all, is the anti-correlatio n between Al and 
Mg (iGratton et al.ll200Tl ICarretta et al.ll2009bl) . It is important to 
mention that such abundance trends are not seen in the halo field 
stars. 

This spectroscopic evidence has been interpreted as the sig¬ 
nature of material processed during H-burning by high tempera¬ 
ture proton-capture reactions (like the Ne-Na and Mg-Al cycles). 
Several theoretical models have been proposed in order to de- 
scribe the formation and early evolution of GCs dP’Ercole et al.l 


I 2 OQ 8 I) . The preferred explanation involves a self-enrichment sce¬ 
nario, within which two subsequent generations of stars co-exist 
in globular clusters and where the second is formed from gas 
polluted by processed material p roduced by massive stars of 
the first (ICaloi & D’ Antonall201 11) . Several sources of processed 
ejecta have been proposed: the slow w inds of intermediate- 
mass AGB_star^D^ntonaeL^SQQ2)j fast ro tating massive 
stars (iDecressin ^ alJ2007^ and massive binaries (Ide Mink et al.l 
20091). A recent att ractive alternative has been proposed by 
Bastian et alJ (1201 3l) . that implies only a single burst of star 


formation. They postulate that the polluted gas concentrates in 
the center and is acquired by low mass stars of the first (and 
only) generation via disk accretion while they are in the pre- 
main sequence fully convective phase. This scenario has the dis¬ 
tinct advantage that it does not require a very large percent¬ 
age of the original cluster population to have been lost. In ad¬ 
dition, variations in heavier elements have also been found in 


Carretta et al 

2010hh. M22 (lMarinoetalH2009). NGC 1851 

Carretta et al 

2011). Terzan5 (Eerraro et alj|2009t). NGC 2419 

Cohen et al. 

2OI0I). M2 (Eardoet all 20131: Yong et all I20l4 

and M75 (jKacharoy. Koch & McWilliam 20131). They are gen- 


erally thought to be the vestige of more massive primitive dwarf 
galaxies that merged with the Galaxy. Therefore, they have im¬ 
portant implicatio ns for the hierar chical merging scenario of 
galaxy formation (iJoo & Leell2013l) . In this paper, we use the 
widely accepted terminology of first and second generation of 
stars referring to the unpolluted and polluted populations regard¬ 
less of the formation scenario. 

NGC 4372 is an old and very metal-poor globular cluster that 
despite its low degree of central concentration, has received little 
attention, mostly due to reddening problems. It is located close 
to the Galactic disk and in a dusty area of the Southern Musca 
Constella tion (RA ^ 1 2:25:45.40, Dec=-72:39:32.4; 1=300.99°, 
b=-9.88° lHarrisl[l996l 2010 edition). At a Galactocentric dis¬ 
tance of Rgc = 7.1 Kpc, NGC 4372 is listed in the Harris cat¬ 
alog as having a low metallicity of [Ee/H] = -2.17 and a fairly 
high reddening of E(B-V) =0.39, which makes it a visually chal¬ 
lenging object. In addition, NGC 4372 h as been claimed to be 
dynam ically associated with NGC 2808 (ICasetti-Dinescu et al.l 

l2^ . 

_ Ea rly photometric studies of this cluster (iHartwick & Hess^ 

II 973 I : lAlcainol 119741 : iBrocato et al.l 1 19961) already reyealed a 
color-magnitude diagram characteristic of yery metal-poor clus¬ 
ters, with a well-defined horizontal branch that extends far to 
the blu e, and with a large and yariable absorption. I Alcaino et al.l 
(I 1991 I) present B VRJ CCD photometry in two oyerlapping fields. 
By comparison with theoretical isochr ones, they deriye an age 
of 15 +4 Gyr. In a more recent study, [^senberg et al.l (l2000l) 
present a homogeneous photometric catalog of 39 GGCs in the 
southern hemisphere. They report for NGC 4372 a foreground 
reddening of E(B-V) = 0.42 , horizontal branch (HB) leyel of 
Whb = 15.30 with a HB ratio (B-R)/(B-hV-hR) = 1.00. They ex¬ 
plain the broadening of the CMD sequences as a consequence of 
the high differential reddening probably due to the nearby Coal- 
sack Nebulae. A de-reddened, narrow er CMD has recently been 
constructed by lKacharoy et alJ (1201 4l) 

Spectroscopic studies of NGC 4372 are limited to 
medium-reso l ution s pectra. Using the near-infrared Can triplet, 
iGeisler et alJ (1 19951) determine a mean metal abundance of 
[Ee/H] = -2.10 + 0.04 from 11 giant stars. They also report 
a mean heliocentric radial yelocity of Vmd - 73.2 +1.4 km/s . 
Using a similar technique, [Rutledge, Hesser & Ste ts^ ( 19971) 
estimate a yalue of [Ee/H] = -2.03 + 0.03 on Zinn & Wesll 
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Fig. 1. A 15'x 15' digitized sky survey image centered on NGC 4372. 
North is up and east is to the left. The red symbols correspond to the 
spatial distribution of the 7 stars analyzed. Th e dashed circles corre¬ 
spond to the core radius, r^= 1.7 ’ (lHarridll99d 2010 edition), and the 
half-light radius, i/,=3.4’ iKacharov et alJ 2014ll . 
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Fig. 2. Color-magnitude diagrams of the globular cluster NGC 4372 on 
an area of 6.0'x 5.2'. The red symbols correspond to the analyzed RGB 
stars'. 


(IT9^ ZW84) scale and a value of [Fe/H]= -1.88 + 0.05 on 
ICarretta & Grattonid 19971 CG97) sc ale. A more up dated value of 
[Fe/H]=-2.19 + 0.08 is provided bv ICarretta et alJ (l2009ah . This 
value is determined from the previously mentioned Ca ii triplet 
data but based on a more accurate metallicity scale. No further 
chemical study has been performed. 

In this paper, we present the first extensive study of NGC 
4372 from high-resolution spectroscopy. We perform an abun¬ 
dance analysis of a large number of elements and analyze the 
stellar population of the cluster. This paper is organized as 
follows: Section 2 describes the observations and data reduction 
and Section 3 describes the methodology used to obtain the 
atmospheric parameters and chemical abundances. Section 4 
presents our results including analysis of iron-peak elements, a- 
eleigents and any (anti-)correlation. Finally, Section 5 presents 
a summary of our main results. 

2. Observations and Data Reduction 

The present work is based on the data collected by the GES 
from the beginning of the Survey up to the end of June 2013. 
These observations are referred to as ^esiDR2iDR 3 and are avail¬ 
able in side the GES collaboration. As explained in iGilmore et al.l 
(12012|) . the GES consortium is structured in several working 
groups, WGs, having specific duties from target selection and 
data reduction to delivery of science data. The data reduction 

' The photometric catalog has been generated by WG5 based on 
UBVI archival images and collected at the Wide-Field Imager (WFI) at 
the 2.2m ESO-MPI telescope (Programes: 164.0-0561 (PI: Krautter), 
68.D-0265 (PI: Ortolan!) and 69.D-0582 (PI: Or tolani)). The images 
were pre-reduced using IRAF package MSCRED ( IValdeslll998h . while 
the stellar photo metry was derived by using the DAOPHOTII a nd ALL- 
STAR programs dStetsonll 19871119921 . Stetson standard fields dStetsonl 
l2000h have been used to photometrically calibrate the data, while as- 
trometric calibration has been p erformed with the c atalogs UCAC3 
dZacharias et alJiTolOh and GSC2 dMcLean et al.ll20(}Qll . 



Teff(K) 

Fig. 3. Parameters of the stars in the Teff- log g plane. The solid line 
corresponds to the isochrone co mputed with the web-tool of the PAR¬ 
SEC group dBressan et al.ll201^ . 


has been performed using a workflow specifica lly developed for 
and by GES that runs the ESQ pub lic pipeline ((Modigliani et al.l 
l2004 iModigliani & Larsenll2012ll . The pipeline has been opti¬ 
mized to reduce GES data and performs automatically sky sub¬ 
traction, barycentric correction and normalization. In addition 
the pipeline calculates radial velocities and a first guess of the 
rotational velocities. A quality control is performed using the 
output parameters from the ESO pipeline, by a visual inspection 
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of the spectra and by the analysis of the signal-to-noise ratio of 
the spectra. A d etailed descr i ption of the data reduction methods 
can be found in ISacco et al.l (l2014l) . 

The present study makes use only of the UVES observa¬ 
tions. While a small number of fibers can be dedicated to UVES, 
their spectra have a higher resolution than GIRAEEE and, just 
as important, have ten times wider spectral coverage. Eor this 
reason, UVES data are crucial for studies of precise chemical 
abundances of a large number of elements. As explained in Sec¬ 
tion [3 this paper makes use only of Concepcion Node atmo¬ 
spheric parameters and chemical abundances rather than the rec¬ 
ommended GES values. Unfortunately, this Node does not ana¬ 
lyze GIRAEEE spectra. Although ideally the used of GES GI¬ 
RAEEE spectra for NGC 4372 would strengthen the conclusion 
of this paper, the different techniques and methods employed, for 
this specific science case, would introduce an undesirable hetero¬ 
geneity. 

In this paper, we present a detailed chemical abundance anal¬ 
ysis of seven red giant branch (RGB) stars observed with UVES. 
The stars are observed in the set-up with central wavelength 580 
nm. The spectra are taken in two arms, resulting in a wavelength 
coverage of 470-684 nm with a gap of ~ 5 nm in the center. The 
ELAMES-UVES fibers have an aperture on the sky of 1”, re¬ 
sulting in a resolving power of R=47000. The observations were 
taken between March 11 and March 16 of 2012 with a mean 
signal-to-noise ratio in the spectrum of ~ 50. Targets for UVES 
were high probability cluster members and in particular evolved 
stars. Eigure[T] shows the spatial distribution of the 7 RGB stars 
covering a wide area around NGC 4372 while Eigure|2]presents 
color-magnitude diagrams of the cluster with the analyzed stars 
indicated. Target selection for globular clusters has been made 
by WG5 (Calibrators and S tandards; coordinator: E. Pan cino). 
More details can be found in IPancino et al.l (12015. in prenll . 

The radial velocity information, available for all the tar¬ 
gets observed, is determined within WG8 by cross-correlation 
against real and synthetic templates. The typical error on the ra¬ 
dial velocities of UVES targets is about 0.4 km s“'. The mean 
radial velocity value for the seven targets is < Vrad > -12.6 ± 1.3 
km/s while the dispersion is 3.6 km/s. The mean r adial velocity 
is in e xcellent agreement with the value found by iGeisler et al.l 
(Il995h of < Vrad >=73.2 +1.4 km/s where the e rror is given by 
the sta ndard error of the mean. More recently, iKacharov et al.l 
(l2014l) present an extensive kinematic study of this cluster using 
ELAMES/GIRAEEE observations of 131 stars. They find a mean 
radial velocity of Vr=75.9 + 0.3 km/s and a central velocity dis¬ 
persion of ctq = 4.5 + 0.3 km/s. As part of t he Gaia-ESO surve y 
the kinematics of seven GGCs is present in lLardo et al.l (l2014l) . 
They obtained from ELAMES/GIRAEEE spectra of more than 
100 stars a mean radial velocity of Vr=75.2 + 0.4 km/s with a 
velocity dispersion of 3.9 km/s, confirming the membership of 
our objects. The stellar parameters of stars analyzed in our study 
are summarized in Table [T] Photometric magnitudes of the ob¬ 
served stars are compiled by GES and presented in Table |2] 

3. Atmospheric Parameters and Abundance 
Analysis 

The astrophysical parameters obtained from the analysis of the 
gesiDR2iDR3 data set will be part of the first Gaia-ESO public 
release of advanced data products. Within WGll, the spectro¬ 
scopic analysis of UVES data is performed by 13 sub-groups, so- 
called Nodes. A multiple parallel analysis of the full data set has 
been implemented where different Nodes use different method¬ 
ologies. Unfortunately, only a small number of Nodes reported 


atmospheric paramete rs for NGC 4372 stars. As explained in 
ISmilianic et al.l (l2014ll . GES recommended parameters are com¬ 
puted only if at least 3 Nodes provided parameters for a given 
star. This decision is made based on internal policy rather than 
the reliability of the values. We decided to make use of Con¬ 
cepcion Node atmospheric parameters and chemical abundances 
during the following analysis. A critical evaluation of the per¬ 
formance of Co ncepcion Node, based on a series of calibrators, 
is presented in (ISmilianic et al.llT014ll . In addition, Concepcion 
Node is the only one that provide, for the 7 stars observed in 
NGC 4372, key elements such as O or Na. 

The atmospheric parameters for NGC 4372 reported by 
the Concepcion Node correspond to the GES iDRl (internal 
data release 1) and were determined by at least 5 different 
Nodes in a system of multiple parallel analysis. The method¬ 
ology and codes used by each Node are d escribed in de¬ 
tail in Appendix A of ISmilianic et al.l (l2014h and they range 
from the classical method of equivalent width (EW) to the use 
of libraries of observed and/or synthetic spectra. This strat¬ 
egy has two main advantages: 1) Ensure that all sources of 
errors are well-understood and quantified, including method- 
dependent effects, 2) Ensure that all types of objects can be 
properly analyzed. To guarantee the homogeneity of the final re¬ 
sults, a number of constraints have been imposed. These con¬ 
straints inclu de: the use of a common line list (iHeiter et al.l 
12015. in nrenl) . the use of one single set of mode l-atmospheres 


(the MARCS models, iGustafsson et ak ( 2008lB. the use of 


single synthetic spectmm lib rary (Ide Lavernv et ^ 12012^ 

S Blanco_et^L[ 2015. in nrenl). a common solar zero point 
;sse. Asplund & SauvaJ2007l) and the analysis of common 
calibration targets. In order to understand the precision and ac¬ 
curacy of the atmospheric parameters, GES makes use of the 
Gaia benchmark stars and a set of calibration clusters. The ac¬ 
curacy is judged by the ability of a given Node to recover the 
reference atmospheric parameters of the benchmark stars from 
their analysis. The precision is judged by the ability of a Node to 
reproduce their own results from multiple analysis of the same 
star. A detail explanation of the performance of each Node as 
well as how the atmospheric pa rameters were determined can be 
found in ISmilianic et al.l (l2014l) . Table[T]summarizes the derived 
atmospheric parameters for NGC 4372 stars. The error reported 
for each parameter corresponds to the dispersion among the re¬ 
sults from different methodologies. This method-to-method dis¬ 
persion is defined as the degree to which multiple methodologies 
can agree on the abundance of a star. Although this dispersion is 
not properly the physical uncertainty of the values, it is a good 
indicator of the precision and is adopted as the typical uncer¬ 
tainty. To further investigate the reliability of these values Eig. |3] 
shows the derived atmospheric parameters of the stars compared 
with those derived from isochrones. 

The chemical abundances for Mg, Al, Si, Ca, Ti, Cr 
and Ni were obtained using EWs of the spectral lines. The 
EWs were determined w ith the automatic code DAOS PEC 
(IStetson & Pancinoll2008l) . GALA (iMucciarelli et al.ll2013h was 
used to determined the elemental abundances. To complete 
the analysis and include other key elements the abundances 
of O, Na, Sc, Y, Ba, and La were also obtained. Eor this 
set of elements whose lines are affected by blending, the 
spectrum-synthesis method were used. Th e local Therm ody- 
namic Equilibrium (LTE) program MOOG (ISnedenlll973l) was 
used for this purpose. Eive synthetic spectra having differ¬ 
ent abundances for each line were calculated, and the best¬ 
fitting value estimated as the one that minimizes the RMS 
scatter. Only lines not contaminated by telluric lines were 
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Table 1. Basic Parameters of the Observed Stars. 



Cname“ 

Vrnrf 

(km/s) 

SNR* 

Teff 

(K) 

log(g) 

(dex) 

(km/s) 

[FelH] 

(dex) 

1 

12250638-7243067 

68.6 ± 0.5 

51.17 

4586 ± 86 

0.70 ± 0.20 

1.56 ±0.10 

-2.23 ± 0.21 

2 

12250660-7239224 

67.9 ± 0.3 

34.63 

4646 ± 56 

1.17 ±0.15 

1.68 ±0.03 

-2.22 ±0.15 

3 

12253419-7235252 

72.7 ± 0.7 

55.73 

4527 ± 33 

0.55 ± 0.07 

1.50 ±0.17 

-2.15 ±0.07 

4 

12253882-7245095 

76.0 ± 0.6 

45.45 

4656 ± 114 

1.19 ± 0.47 

1.56 ± 0.20 

-2.21 ± 0.20 

5 

12260765-7236514 

75.2 ± 0.5 

71.14 

4466 ±38 

0.55 ± 0.09 

1.65 ±0.05 

-2.20 ±0.17 

6 

12264293-7241576 

71.0 ±0.7 

44.35 

4712 ± 87 

1.19 ±0.21 

1.55 ±0.13 

-2.15 ±0.12 

7 

12264875-7239413 

76.9 ± 0.4 

43.89 

4633 ± 68 

0.84 ± 0.34 

1.54 ±0.05 

-2.44 ±0.16 


“ GES object name. It is formed from the coordinates of the object splicing the RA in hours, minutes and 
seconds (to 2 decimal places) and the Dec in degrees, minutes and seconds (to 1 decimal place) together, 
including a sign for the declination. 

* Mean signal-to-noise ratio in the spectrum. 


Table 2. Photometric Properties of the Observed Stars. 



Cname“ 

B* 

(mag) 

k* 

(mag) 

(mag) 

7* 

(mag) 

/A 

(mag) 

/C 

(mag) 

E{B - vy 

(mag) 

(mas/yr) 

Boec 

(mas/yr) 

1 

12250638-7243067 

14.70 

13.22 

11.52 

10.30 

9.61 

9.44 

0.57 

-6.10 

4.80 

2 

12250660-7239224 

15.22 

13.80 

12.11 

10.86 

10.15 

9.98 

0.57 

-12.40 

6.70 

3 

12253419-7235252 

14.83 

13.21 

11.29 

9.92 

9.18 

8.97 

0.59 

-7.10 

8.60 

4 

12253882-7245095 

15.15 

13.84 

12.21 

11.10 

10.43 

10.28 

0.54 

-5.50 

9.90 

5 

12260765-7236514 

14.66 

13.12 

11.25 

9.92 

9.19 

8.97 

0.58 

-13.80 

9.00 

6 

12264293-7241576 

15.08 

13.76 

12.14 

10.98 

10.30 

10.16 

0.54 

-1.10 

0.60 

7 

12264875-7239413 

15.16 

13.82 

12.08 

10.87 

10.17 

9.98 

0.55 

-0.10 

13.40 


“ GES object name. It is formed from the coordinates of the object splicing the RA in hours, minutes and seconds (to 
2 decimal places) and the Dec in degrees, minutes and seconds (to 1 decimal place) together, including a sign for the 
declination. 

Magnitudes from APASS ( AAVSO photometric All-sky Survey). 

Magnitudes from 2MASS iSkrutskie et alJ[20()3l . 

Reddening values from the Galactic dust extinctionjSchlegel ^Ejnkbeiner & Davislll99^ . 

'■ Proper motion in RA and Dec from UCAC dZacharias et alj2010h . 


Table 3. Abundance Ratios. 


El. 

Method" 

Nlines* 

1 

2 

3 

4 

5 

6 

7 

Mean" 

[0/Fe] 

SYNTH 

1 

0.22 

0.36 

0.53 

0.46 

0.67 

0.39 

0.80 

0.44 ± 0.07 

[Na/Fe]MLTE 

SYNTH 

2 

0.39 

0.91 

0.53 

0.69 

0.07 

0.67 

0.53 

0.54 ±0.13 

[Mg/Fe] 

EW 

4 

0.41 

0.69 

0.37 

0.45 

0.20 

0.53 

0.62 

0.44 ± 0.07 

[Al/Fe] 

EW 

2 

0.65 

1.16 

0.85 

1.12 


1.01 

0.01 

0.96 ±0.10 

[Si/Fe] 

EW 

1 



0.46 

0.52 

0.48 


0.77 

0.48 ± 0.02 

[Ca/Fe] 

EW 

18 

0.25 

0.29 

0.19 

0.24 

0.22 

0.28 

0.60 

0.24 ± 0.02 

[Sc/Fe] 

SYNTH 

1 

-0.04 

-0.09 

-0.14 

-0.14 

-0.03 

-0.11 

-0.01 

-0.09 ± 0.02 

[Ti/FeY 

EW 

14 

0.29 

0.40 

0.24 

0.30 

0.23 

0.38 

0.53 

0.31 ±0.03 

[Fe/H] 

EW 

100 

-2.23 

-2.22 

-2.15 

-2.21 

-2.20 

-2.15 

-2.44 

-2.19 ±0.02 

[Cr/Fe] 

EW 

6 

-0.38 

-0.20 

-0.33 

-0.37 

-0.36 

-0.28 

-0.15 

-0.32 ± 0.03 

[Ni/Fe] 

EW 

6 

0.07 

-0.03 

0.14 

-0.03 

-0.06 

0.07 

0.06 

0.03 ± 0.03 

[Y/Fe] 

SYNTH 

1 

0.05 

0.01 

0.03 

0.00 

0.10 

0.09 

0.06 

0.05 ± 0.02 

[Ba/Fe] 

SYNTH 

1 

-0.29 

-0.15 

-0.09 

-0.16 

-0.11 

-0.17 

-0.23 

-0.16 ±0.03 

[La/Fe] 

SYNTH 

1 

0.09 

-0.03 

-0.36 

-0.31 

-0.06 

-0.06 

0.14 

-0.12 ±0.08 


“ Method used in the abundance determination: SYNTH (Spectrum-synthesis method), EW (Equivalent 
Width method). 

* Average number of lines employed in the abundance determination. 

Values obtained excluding star #7. 

Values obtained from the average of Ti i and Ti ii. 


used. Eor both methods, the GES guidelines rega rding the 
use of the GES line list (iHeiter et al. 2015. in prepi). th e use 


of MARC S model atmospheres (iGustafsson et al. 


2008h and 


the use of iGrevesse, Asplund & Sauvall ( 200% solar-zero point 
were strictly followed. In addition, the Na abundance was re¬ 
computed and a non-LTE correction was applied based on 


iMashonkina, Sbimanskii & SakhibullinI (l2000h . All the used 
Na abundances are N LTE corrected. We refer the reader to 
ISmilianic et al.l (1201 4|) for a detailed description of the abun¬ 
dance analysis. The chemical abundances of the observed stars 
are presented in Table [3] 
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Fig. 4. Abundance ratios versus [Fe/H] for individual members of the cluster. The dashed lines indicate the Icr area around the average value. 


An internal error (star-to-star) analysis of the derived abun¬ 
dances has been performed to properly quantify the internal 
spread of abundance within a cluster. 

Table |4] lists the two primary sources of errors contribut¬ 
ing to the total budget {(Trot)- the uncertainties in the measure¬ 
ments and the uncertainties associated with the atmospheric pa¬ 
rameters. The parameter ctew corresponds with the mean chem¬ 
ical abundance variation due to error on the EW measurements. 
Uncertaintie s in the measurements o f EWs are computed by 
DAOSPEC dStetson & Pancino||2008l) . These confidence inter¬ 
vals estimate the goodness of the fit but also take into account 
the quality of the spectrum (resolution, S/N, s pectra defects,...). 
These EW errors , EWerr, are used by GALA dMucciarelli et alJ 
l2013h to obtain the uncertainty on the abundance of each element 
by varying the EW by I EWerr- For the elements whose abun¬ 
dance was obtained by spectrum-synthesis ctew was obtained as 
the error given by the fitting procedure. An internal error analysis 


was performed by varying Tgff, log(g) and ^ and redetermined 
abundances of star #3, selected as having representative atmo¬ 
spheric parameters. Parameters were varied by ATe/f — h- 70 K, 
Alog(g)= -1-0.19 and A^=-i-0.14 km s“*, which correspond with 
the mean error of the sample parameters. We note that these 
parameter variations are smaller than the standard steps in the 
MARCS model grid. The MARCS model interpolator^ has been 
used to generate intermediate models from the initial grid. The 
program interpolates the thermal structure (T), the electronic 
pressure (Pg), the gas pressure (P^), the opacity {k) and the micro 
turbulence velocity (^) as a function of Tgff, log(g) and metallic- 
ity [M/H]. The interpolator has been extensively tested on a pre¬ 
vious grid of MARCS models considering the following range 
of parameters: 3800 K < Teff < 7000 K, 0.0 < log g < 5.0, -4.0 < 


^ The program and a detailed user manual are available on the MARCS 
web site: http://marcs.astro.uu.se 
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Table 4. Sensitivity of Derived Abundances to the Atmospheric Parameters and EWs. 



ATeff^lQ (k) 

Alog g=0.19(dex) 

A^=0.14 (km/s) 

m EW 

^ tOT 

^obs 

A[6)/Ee] 

0.04 

0.07 

0.03 

0.04 

0.09 

0.15 

A[Na/Fe] 

0.00 

-0.07 

-0.03 

0.04 

0.09 

0.29 

A[Mg/Fe] 

0.11 

-0.08 

-0.07 

0.02 

0.15 

0.16 

A[Al/Fe] 

0.02 

0.04 

0.03 

0.04 

0.07 

0.21 

A[Si/Fe] 

0.04 

0.00 

-0.02 

0.02 

0.05 

0.03 

A[Ca/Fe] 

0.08 

-0.01 

-0.06 

0.02 

0.10 

0.04 

A[SclFe] 

0.12 

0.16 

0.10 

0.05 

0.23 

0.05 

A[TilFe] 

0.00 

0.07 

-0.04 

0.05 

0.09 

0.07 

A[CrlFe\ 

0.14 

0.01 

-0.10 

0.03 

0.17 

0.07 

A[NilFe] 

0.13 

0.01 

-0.08 

0.02 

0.15 

0.08 

A[YIFe] 

0.03 

0.07 

-0.01 

0.06 

0.10 

0.04 

A[BalFe\ 

-0.09 

-0.06 

-0.02 

0.06 

0.13 

0.07 

A[LalFe\ 

0.12 

0.14 

0.08 

0.07 

0.21 

0.17 


Note. - The sensitivity determination was performed for star #3. 

(Tew is the error in the measurements. 

(Jtot is the squared root of the sum of the squares of the individual errors. 
(Tobs is the mean observed dispersion. 



Fig. 5. Box plot of the NGC 4372 star element abundances. Star #7 has 
been excluded. For each box the central horizontal line is the median 
of the data, while the lower and upper lines indicate the first and third 
quartile. The whiskers extend out to the maximum and minimum value 
of the data. Outliers are identified with small circles. 


[M/H] < 0.0 (lMasseronll2006h . The interpolation is optimized to 
account for non-linearities in the grid and the new interpolated 
model must lie inside a complete cube of existing models in the 
parameter space (Tg//, log g, [M/H]). With the actual grid pa¬ 
rameter steps, maximum errors in the interpolated quantities re¬ 
main below 0.25% and a few % for Pg and Pg. Within these lim¬ 
its, we do not consider that the error budget might increase dra¬ 
matically and that the mean method-to-method dispersion drives 
the sensitivity to atmospheric parameter variations. This estima¬ 
tion o f the internal error w as performed following the prescrip¬ 
tion of iMarino et al.l (l2008h . The final total error, (cr,g„), has been 
computed as the square root of the sum of the squares of the 
individual errors. Table|4]lists also the observed star-to-star dis¬ 
persion {(Tohs)- We remark that our goal is to search for evidence 
of star-to-star intrinsic abundance variation in each element by 
comparing the observed dispersion and the internal errors 

(cTtot)- For this reason, external sources of error as systematics 
that do not affect relative abundances are not considered. 


4. Results 

4.1. Chemical homogeneity of the cluster 

Chemically inhomogeneous populations are observed in vir¬ 
tually all massive, old globular clusters well-studied to date, 
but no t in open clusters, with the possible except i on of NGC 
6791 (iG eisler et ap 1201^ but see [Bragaglia et al.l (l2014ll and 
ICunha et al.l (1201 41 ) 1 . Studying the chemical homogeneity or in¬ 
homogeneity of a cluster is necessary to better understand the 
mechanism of their formation. A large body of evidence now 
shows that t he stars in a GC do not share the s ame chemical 
composition (iGratton, Carretta & Bragagliai2012h . As a general 
rule we can define a globular cluster as an object homogeneous 
in its Fe content and most other heavy elements, but the light 
elements Li, C, N, O, Na, Mg and Al can show substantial intr¬ 
acluster variations. We investigate the degree of inhomogeneity 
of the abundances of cluster members in NGC 4372. 

Figure 0 ] shows the abundance ratios versus [Fe/H] for indi¬ 
vidual member stars in the cluster. In each panel, the intersec¬ 
tion of the dashed lines delimit a Icr area around the average 
value. Star-to-star error bar is indicated. Star #7 stands out as 
slightly more metal poor than the main body of cluster mem¬ 
bers. The a-elements of this peculiar star also stand out from 
the cluster mean abundances, except for Mg. The radial veloc¬ 
ity and stellar parameters for this object shown in Table [T] are 
in agreement with those of the rest of the observed stars but the 
distinct chemical pattern of this star suggests wither that it has 
some kind of anomaly or it may not be a member. The possi¬ 
bility of star #7 being in a binary system would likely alter its 
radial velocity as well as artificially lower its derived metallicity 
due to increased continuum flux. As part of the quality control, 
GES final products include two binary flags: a) a visual inspec¬ 
tion of the cross-correlation function (CCF), computed before 
co-adding multi-epoch observations, has been performed. A star 
is flagged as a candidate double-lined spectroscopic binary, if the 
CCFs are characterized by the presence of more than one peak 
or a single peak with strong asymmetries; b) a star is classified as 
a single-lined spectroscopic binary if the median absolute devia¬ 
tion of multi-epoch repeated measurements of the radial velocity 
is larger than twice the error on the radial velocity. None of the 
NGC 4372 observed stars have been flagged as binary candi¬ 
dates. Given its spatial location, velocity, position in the CMD 
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and low metallicity^ it see ms very unlikely that this is a field 
star. iKacharov et al. (l2014t) show that according to the predic¬ 
tion of the Besancon Galactic model, only a few field stars with 
that velocity are expected in the direction of NGC 4372 where 
none of them are more metal-poor tha n -1.8 dex. This suppo rt the 
view that star # 7 is a cluster member. lLapenna et alJ ( 1201 4|) sug¬ 
gest that NLTE effects driven by over ionization mechanisms are 
present in the atmosphere of AGB stars and significantly affect 
Fel lines, while leaving Fell features unaltered. This effect could 
underestimate the metallicity of AGB stars. The low metallicity 
value of this peculiar star could be due to the fact that it is an 
AGB star. It is still unclear why the a-elements of this star are 
also significantly different. To be cautious and conservative, we 
exclude star #7 from the analysis. 

Fig Elpresents the abundance pattern analyzed in NGC 4372. 
The box plot illustrates the median and the interquartile range 
(IQR) of the derived values. Possible outliers are also included 
where an outlier is defined if it deviates by more than 1.5 IQR. 
A large abundance ratio range is present for O, Na, Mg and Al. 
The star-to-star variations are smaller for the heavier elements. 
Fa presents a large abundance range but we note that the errors 
in the derived abundances of Sc and Fa are very large. For both 
elements, the spectrum-synthesis method with only a single line 
(Sen /15684.202 A and Fan 45122.995 A) was used to derived 
the abundances. The sensitivity of the derived abundances for 
these two elements to the atmospheric parameters, and specif¬ 
ically to Teff and log g (see Table IH, is very significant and 
any conclusion drawn from these two elements should be treated 
with caution.To further investigate the degree of homogeneity 
of the abundances, we compare the scatter produced by inter¬ 
nal errors, (cr,g,), with the observed dispersion in the chemical 
abundances, {(Tobs)- These values correspond to columns 6 and 7 
of Table |4l One can consider inhomogeneity when the intrinsic 
scatter is significantly higher than the expected dispersion given 
by the internal errors. If we exclude star #7 from the error esti¬ 
mation, a clear intrinsic spread can be identified only in O, Na 
and Al. 



-3 -2 -1 0 

[Fe/H] 


Fig. 6. Iron-peak element ratios ([Sc/Fe], [Cr/Fe], [Ni/Fe]) ver¬ 
sus [Fe/H]. Red symbols correspond to values for NGC 4372 
from the present study where the open red circle corresponds to 
the peculiar star #7. Blue stars are Galactic Globular Clusters 


divans et alJl2001l:lLee. Camev & Habeoodll2005l:lLee & Camevll2002: 

Carrettal 20061; ICarretta et aLll2009al 

2010al: IVillanova & Geisleill201 1 

Munoz. Geisler & Villanoval 1201.31: 

Koch & MeWi 

liami 20111). Grav 

symbols are halo and disk stars dFulbrisht 20001; 

Reddy et all 120031; 

Reddv. Lambert & Allende Prietol2006l:ICavrel et a 

.|2004|) while black 


symbols correspond to extra-galactic obje cts: Draco, Sextans, Ursa 
minor, Carina and Sagittar i us dwarf galaxy dShetrone. Cote & Sargent! 
l2001l : ISbordone et al.l200'^ : lKoch et al.l20()8 al) and th e ultra-faint dwarf 
sphero idals Bootes I and Hercules dlshisaki et all l2014l : iKoch et ^ 
l2008bh . Star-to-star error bar is indicated. 


4.2. Iron and Iron-peak elements 

We obtain a mean metallicity for NGC 4372 of [Fe/H]= -2.23 + 
0.10. Figure. |4] shows no evidence for an intrinsic Fe abundance 
spread with the exception of star #7. As discussed in the previ¬ 
ous section, this star is anomalous in its chemical behavior and 
we prefer to exclude it from cluster means. The issue of intrin¬ 
sic metallicity spreads in GGCs is of great current interest. Such 
spreads are found generally only in the most luminous GGCs, 
with My < -8.5. NGC 4372 has My = -7.8 so it is unlikely to 
host an intrinsic Fe abundance spread. However, we do note that 
iGeisler et alJ d 19951) did find one of their sample of 11 stars to be 
significantly (about 0.5 dex) more metal-poor than the rest. This 
star, like star #7, also had a velocity compatible with member¬ 
ship. Excluding star #7, we found a mean [Fe/H] value of: 

[Fe/H]=-2.19 ±0.03 

The first attempt t o deriv e a metallicity for NGC 4372 
was by IZinn & West! (11984 ZW84) o btaining a v alue of 
[Fe/H]= -2.08 + 0.15. Most recently, IGeisler et al.l d 19951) 
analyzed medium-resolution spectra of ten giant stars and 
obtain, through the near-infrared Ca ii triplet technique, a 
mean metallicity of [Fe/H]=-2.10 + 0.04. Several other au¬ 
thors have attempted to derive the metallicity based on 
similar techniques but using different scales with a variety 


of res ul ts (eg. [Fe/H1zvyR4=-2 .03; [Rutledge, Hesser & Stetson! 
dl997l) '). ICarretta et al.l d2009cl) adopt a new scale that is a fun¬ 
damental shift from the older and widely used ZW84 metallicity 
scale. The authors argue that this traditional scale was calibrated 
against only a handful of hi gh-dispersion spectro scopic [Fe/H] 
values available at that time. Carretta et al.l d2009cl) define an ac¬ 
curate and updated metallicity scale using high-dispersion and 
high signal-to-noise spectroscopic measures of 19 GCs cover¬ 
ing the metallicity range of the bulk of GGCs. Based on this 
scale they provide a value of [Fe/H ] =-2.1 9 ± 0.08 for NGC 4372 
metallicity, based on IGeisler et alJ (1 19951) Can triplet data. This 
updated result is in excellent agreement with our metallicity. 
The RMS scatter of our met allicity is a measure o f the intrinsic 
spread of iron in the cluster. ICarretta et al.l (l2009cl) confirm that 
the scatter in Fe of most GCs is very small with an upper limit 
of less than 0.05 dex. Our observed scatter is consistent with that 
expected from errors and thus we conclude that there is no clear 
metallicity spread in NGC 4372, with the caveat that star #7 is 
considered either an outlier for other reasons or a non-member. 

The chemical abundances for the iron-peak elements Sc, 
Cr and Ni are listed in Table [3 The values are solar within 
the errors except for the Cr abundance which is underabun- 
dant. Figure |6] shows the elemental abundance of each star 
compared with a variety of galactic and extra-galactic ob¬ 


jects. We have included values from GGCs (Ivans et aT[ 2001 
(Lee. Carney & HabgoodI 120051 : iFee & Carnev 2002; Carrettal 
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20061 Carretta et al J l2009a 

2010al Villanova & Geisle 

20111; 

Munoz. Geisler & Villanoval 

20131 

iKoch & Me William 

20111); 

disk and halo stars (Fulbrightj 

120001; Reddv etal. 

120031: 

Reddv, Lambert & Allende Prietoil2006l: Cavrel et alJ 2004) and 


extxa-galactic objects such as Draco, Sextans, Ursa Mi 
nor and Sagittarius dwarf galaxy and the ultra-faint dwarf 
spher o idals Bootes I and H e rcules d^e trone. C ote & Sargent 
20011 Sbordone et alJ l2007t llshigaki et alJ I 2 OI 4 I: iKoch et al ' 
2008 ^ 


In general, we found that NGC 4372 stars have abundances 
of these elements which agree with those of other GCs and halo 
field stars of similar metallicity. 


4.3. a elements 

All the a elements listed in Table [2 (Mg, Si, Ca, Tij^ are 
overabundant relative to the Sun. This is a common feature 
among almost every GC as well as among similarly metal- 
poor halo_2pld_st^in_the Galaxy. A glaring exception is Rup 
106 (IVillanova et al.llT013h . which shows solar a element abun¬ 
dances. Figure I 2 ] shows the a-element (Mg, Si, Ca, Ti) ratios as 
a function of metallicity. For comparison purpose, we have in¬ 
cluded values from GGCs, disk and halo stars and extra-galactic 
objects. The sources of the data are the same as those given in 
the previous section. The a elements in NGC 4372 seem to fol¬ 
low the same trend as GGCs and are fully compatible with halo 
field stars. Star #7 stands out from the cluster behavior, suggest¬ 
ing once again the singularity of this star. Excluding star #7, we 
derive for NGC 4372 a mean a element abundance of: 

[ajFe] = h-0.37 + 0.07 

Figure[2represents the [a/Fe] versus [Fe/FI] relation. Differ¬ 
ent symbols and colors are defined as in Figure|7] Similar to halo 
Milky Way stars, GGCs show a constant overabundance of a ele¬ 
ments over a wide range of metallicities ([Fe/H] < -1). Contrarily 
metal-rich dwarf spheroidal galaxies, [Fe/H] > -2, tend to have a 
much lower a content than Galactic objects at si milar metallic- 
ity le ading in some cases to even sub-solar ratios (iGeisler et al.l 
l2007h . NGC 4372 falls in a region where both Galactic and ex- 
traGalactic objects overlap in their a element content so it is not 
possible to draw conclusions regarding its origin from this dia¬ 
gram. 


4.4. Na-0 anticorrelation 

The Na-O anticorrelation is the classical signature of the pro¬ 
cess of proton-capture reactions in H-burning at high tempera¬ 
ture in a previous generation of stars. All GGCs studied in detail 
to da te show this feature w ith the possible exception of Ruprecht 
106 (IVillanova et al.ll^l3h . It is important to mention that such 
abundance trends are not seen in the halo field stars with similar 
metallicity. As described in the Introduct ion, this chemical sig - 
nature has been proposed to define a GC (ICarretta et al.ll2010a) . 
Figure |2 shows the Na-O abundances in the stars of our sam¬ 
ple. For comparison purposes, abundances of GGCs and halo 
and disk field stars have been overplotted. In addition, two dis¬ 
tinct Na-O anticorrelations have been over-plotted as blue solid 
lines, which correspo nd to the dilution models determined by 
ICarretta et al.l (l2009bl) . One is O richer and is represented by the 
trend of the stars in NGC 7078 (M15) ([Fe/H] =-2.31) while the 

^ Since O shows a star-to-star variation and the Na-O anti-correlation 
(see Section l4!4l l. it will be treated separately. 



[Fe/H] 


Fig. 7. u-element ratios ([Mg/Fe], [Si/Fe], [Ca/Fe], [Ti/Fe]) versus 
[Fe/H]. Symbols are as in Fig.j^ 
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Fig. 8. Mean a-element abundance versus [Fe/H]. Symbols are as in 

Fig.m 
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other is O poorer and is represented by the trend of the star s in 
NGC 2808 ([Fe/H]=-1.51). Following l^etta et al.l (l2009bll we 
have obtained the dilution model for NGC 4372. The red dash- 
dotted line was obtained using the full sample and the dashed 
line excludes star #7. NGC 4372 stars have a clear intrinsic 
dispersion and apparent anti-correlation, which corroborate the 
chemical inhomogeneities found in section 4.1. Figure |9] con¬ 
firms the initial assumption that NGC 4372 is a multiple popu¬ 
lation GC. Star #1 poses a value of [O/Fe] significantly smaller 
than the rest of the sample. This star does not follow the Na-O 
anti correlation describe by the rest of the stars in NGC 4372. 
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Fig. 9. [0/Fe] vs [Na/Fe]. Red symbols correspond with this study 
where the open red symbol corresponds t o the peculiar star #7. Black 
symbols correspond with GGC data from ICarretta et al] ( l2009ah while 
gray s y mbols correspond w i th halo and disk field stars from | Fu lbrighll 
l2000h : iReddv et ^ (120031) . I Reddy, Lambert & Allende P rietoUgOO^ 
Blue so lid lines correspond to the dilution model from ICarretta et al.l 
( l2009bt) for NGC 7078 (M15) and NGC 2808. The red dash-dotted line 
is the dilution model for NGC 4372 considering the full sample and the 
dashed line excluded star #7. The black lines are the empirical separa- 
ti ons into primo r dial, in termediate and extreme populations according 
to lCarretta et alj ( l2009bt) . Star-to-star error bar is indicated. 


The O abundances were obtained using the line /16300.30 A. 
This line is significantly weaker in star #1 than in the rest of 
the sample. In addition, NGC 4372 stars inhabit an area in the 
figure that follows the general GCs trend, although our sample 
lies at the high [O/Fe] end. 

Some interesting characteristics can be identified among our 
sample. One star is very Na-poor/O-rich, which corresponds to 
the putative primordial stellar component, while a group of Na- 
rich/O-poor stars would be associated with a second generation 
of stars. However the low number statistics can not confirm that 
split. Our sample does not seem to show any star with very large 
O-depletion which would lie in the Extreme region, resembling 
the behavior of other metal-poo r clusters such as NG C 7078 
(M15), NGC 7099 or NGC 4590 dCarretta et alJl2009hh . 

The dilution model considered by ICarretta et al.l (l2009bh to 
explain the Na-O anti-correlation (Fig. |3 blue and red lines), 
makes the basic assumptions: a) the polluting material has a 
well-defined composition and b) that polluting material is then 
diluted with a variable amount of primordial material producing 
the characteristic pattern of the Na-O anti-correlation. The min¬ 
imum Na and maximum O abundances in each cluster represent 
the original Na and O composition of the cluster. NGC 4372 
ha s abundances f or O sl ightly high compared with the sample 
of ICarretta et al.l (l2009bl) . High values of O content (as an a- 
element) would imply a marginal contribution by type la SNe 
to their original composition. However, the rest of the a abun¬ 
dances (Mg, Si, Ca and Ti) analyzed in NGC 4372 are full com¬ 
patible with GGCs (Fig|7]and[8]l. The apparent offset of the Na-O 
anti-correlation observed in Fig.^could indicate that NGC 4372 
was formed in an environment with high O for its metallicity. 


Fig. 10. [Mg/Fe] vs [Al/Fe]. Symbols are as in Fig.|3 Star-to-star error 
bar is indicated. 


*☆ 

• .A 


: -'is#- : 

t - T_' . »• 

. a= 0.87 ± 0.18 

. a= 1.67 ± 0.29 

r =0.94 . . ★ 

r = 0.94 . , . ★ 

•*•?;**■ 

' . * •' 5 . 






T 

.. . . 1—J—1 


- 0.5 0.0 0.5 1.0 - 0.5 0.0 0.5 1.0 

[Al/Fe] [Mg/Fe] 


Fig. 11. Same as Figure[T0]but for [Si/Fe] versus [Al/Fe], [Si/Fe] versus 
[Mg/Fe], [Na/Fe] versus [Al/Fe] and [Na/Fe] versus [Mg/Fe]. Symbols 
are as in Fig.|9] Star-to-star errors bars are indicated in each panel. Dot¬ 
ted red lines correspond to the linear fit excluding star #7. The slope 
of each fit is represented as a with the correspondent sigma error. The 
parameter r corresponds to the Pearson correlation coefficient. 


4.5. Mg-AI cycle 

In addition to the CNO and Ne-Na cycle, there is evidence that 
the Mg-Al cycle is also active in GC polluters. Models predict 
that Al should show correlations with elements that are enhanced 
by the action of the Ne-Na (such as Na) and Mg-Al cycles 
and should anti-correlate with elements that are depleted in H- 
bu rning at high temperature (such as O and Mg) (see discussion 
in iGratton, Carretta & Bragaglia|l2()T^ . We analyze the abun¬ 
dances of some proton-capture elements to identify any (anti- 
)correlation in the cluster. 

The pattern of abundances of the elements participating in 
proton-capture reactions (Na, Al, Mg, Si) observed in NGC 4372 
giant stars is summarized in Figure [10] and Figure [TT] The level 
of scatter in all these plots is significant. For comparison pur¬ 
poses, abundances of GGCs and halo and disk field stars have 
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been over-plotted. If the temperature is high enough, one should 
expect some degree of anti correlation between Mg and Al. How¬ 
ever the production of Al at the expense of Mg via the MgAl- 
cycle does no t result into such a well-d efined anticorrelation as 
the Na-O one (iLanger & Hoffmanll995h . In fact, large variations 
in Al are of ten accompanied by much smaller changes in Mg 
abundances dCarretta et alJl2009a) . The Mg-Al behavior in NGC 
4372 (Figure [TOll prese nts similar features to the ones found in 
studies of other GGCs dCarretta et al.ll20d^ . 

A clear Al spread is present in our sample with a star-to-star 
variation range in the Al abundances of ~ 0.5 dex (excluding 
star #7). However we need to c onsider the potential in fluence of 
NLTE effect in Al abundances. iThvgesen et al] d2014t) show that 
there is a non negligible NLTE effect on at least some Al lines. 
This NLTE effect correlates with Tg// where cooler stars would 
show stronger negative NLTE correction that warmer ones. Tak¬ 
ing into account the small range of stellar parameters (and in 
particular Tg//^) of our sample stars, the Al spread should not be 
affected significantly by NLTE corrections. On the other hand, 
the larger uncertainty on Mg suggests that Mg may be consis¬ 
tent with no spread. The correlation detected among stars of 
NGC 4372 follows the behavior displayed by GGC stars with 
a clear distinction from disk and halo field stars. Unlike the 
well-known O-Na anti correlation, the Mg- Al anti correlation 
is more difficult to reprod u ce in simulations dPenissenkov et al] 
Il998t[^ntura et alJl20Qlh . lCarretta et al.l (l^09a ) show that Al- 
rich and Mg-depleted stars are present only in massive clusters 
(NGC 2808, NGC 6388, NGC 6441), metal-poor clusters (NGC 
6752) or both (NGC 7078 = M15). In those clusters a clear Mg- 
A1 anti correlation is o bserved even among main sequence stars 
dBragaglia et alJl2010h . More precise Mg values would be re¬ 
quired to clarify the present of a Mg spread in NGC4372. The 
presence of a Mg-Al anti correlation in NGC4372 would indi¬ 
cate that the polluted generation has been enriched by material 
from stars where the Mg-Al burning cycle was active. 

Top panels in Eigure[TT]show the abundances of elements in¬ 
volved in the Ne-Na, Mg-Al cycles of proton-capture reactions 
in high temperature H-burning. The dotted red line in each panel 
of Eigure [TT] corresponds to a linear fit of the sample (exclud¬ 
ing star #7) where the parameter a represents the slope of the 
fit with the sigma error. T he Pearson correlation coefficient, r, is 
also shown in each panel. I^rretta et akl (l2009ah found that stars 
with extreme Al overabundance also show Si enhancement with 
respect to the remaining stars in the same cluster. Again, this 
effect is limited to massive or metal-poor GGCs. The correla¬ 
tion between Si and Al abunda nces is a signature of production 
of ^*Si from the Mg-Al cycle (lYong et alll2005h . The reaction 
producing this isotope becomes predominant in the Mg-Al cy- 
cle when the temperature is very h igh (exceeding Tg = 65 K) 
dArnould, Gorielv & Jorissenlll999h . The chemical pattern ob¬ 
served in these giant stars must be imprinted by a previous gener¬ 
ation of massive stars to be capable to reach those high tempera¬ 
tures. If we exclude star #7 from the analysis (open red symbol). 
Si does not show large star-to-star variations but unfortunately 
the number of stars with reliable Si abundances are very limited. 
A clear Na-Al correlation is evident in our data. This is a reflec¬ 
tion of the relation between Mg-Al and Ne-Na cycles (see dis¬ 
cussion in lGratton, Carretta & Bragagliall20f^ . In fact, Al and 
Na are predicted to be simultaneously enhanced when the Ne- 
Na and Mg-Al cycles are both acting. A Na-Mg relation is also 
identified in NGC 4372 stars. This behavior does not follow the 
general trend of GGCs and in fact one would expect to find a 
Na-Mg anti correlation. This relation is driven mostly by two 
stars only. Large Mg errors could cause this unexpected relation. 



Atomic Number (Z) 

Fig. 12. Abundance distribution of heavy neutron-capture elements in 
NGC 4372. Black symbols correspond to the individual stars while red 
symbols represent the mean abundance for the cluster (excluding star 
#7). Solid and das hed lines show the sca led-solar r- and ^-process pat¬ 
terns respectively dSimmerer et alll2004l) and have been normalized to 
the mean abundance of Ba. Notation: log e (El)= logio(N£//Na) + 12.0 


More extensive samples of stars with accurate determinations of 
abundances are required for a more detailed analysis. 

4.6. Heavy Elements 

While light-element variations are well-known in GCs, intrinsic 
dispersion among heavier elements is less common. Most of the 
heavier elements (Z > 30) are produced either by slow or rapid 
neutron-capture reactions (the so-called s- and r- processes). S- 
process happens in a different physical condition with respect to 
r-process and are thus likely to happen in different astrophysical 
sites. We measured the abundances of the neutron-capture ele¬ 
ments: Y, La and Ba. These elements are mainly produced by the 
i-process at solar metallicity in which the neutron-capture time 
is much longer than the beta-decay lifetime. Erom our data we 
could not determine the r- and s-process ratio of neutron capture, 
because we did not measure a reliable abundance for a typical r- 
process species (e.g. Eu). 

The abundance pattern of heavy neutron-capture elements 
(Y, Ba, La) is shown in Eig. fl^ We have over plotted the pure 
r- and s-process patterns from Simmerer et alJ (l2004l) . The Ba 
abundance was used to set the zero point of the curves. Although 
the analysis includes very few points for any solid conclusion, 
the observed abundance of Y can be reproduced by a pure r- 
process pattern. 

Eigure [13] shows the abundance ratios of the heavy elements 
analyzed as a function of [Na/Ee] and [Al/Ee]. The [Ba/Ee] 
abundance ratio is sub-solar in agreement with recent studies 
which found that [Ba/Ee] decreases with low metallicites and 
beco mes sub-solar at least in fields stars wi th [Ee/H] < -2.0 
(e.g. iGratton & Snedenlll994j: lEulbrigh3l2002h . These sub-solar 
values will agree with recent studies that indicate that metal- 
poor globular clusters primarily exhibit r-process signatures, 
due to the inefficiency of ^-process in low metallicity environ¬ 
ments d Gratton, Sneden & Carr ettal [20041 [Roederer et al.ll20Tol 
[Koch, Cote & McWilliamll^09ir At [Ee/H] < -2 dex, the envi¬ 
ronment was not yet polluted by a sufficient number of low mass 
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Fig. 13. Neutron capture elements [Y/Fe], [Ba/Fe] and [La/Fe] as a function of light elements [Na/Fe] and [Al/Fe]. Open symbol corresponds to 
star #7. 



Fig. 14. Relation of the light and heavy ^-process elements, [hs/ls\ = \hdi[Y}, as a function of [Ba/Fe]. Open symbol corresponds to star #7. The 
dashed line corresponds to the linear fit excluding star #7. The solid line corresponds with the 1:1 correlation. The slope of the fit is represented as 
a with the correspondent sigma error. The Pearson correlation coefficient is also shown. 


AGB stars and the average s-process abundances are low. All the 
elements show an apparent small dispersion however, the ob¬ 
served range is very comparable to that expected from the errors 
(see Table so we can not confirm any clear evidence for an 
intrinsic spread in any of these elements. 


Only a handful of metal-poor globular clusters show 
a pote ntial star-to-sta r dispersion in neutron-capture ele - 
ments ( RoedereJ 20111: iKacharov. Koch & McWilliamI 1201^ . 
iMarino et al.l ( 2009 ) found a wide range of abundances val¬ 


ues for s-process elements, Y, Zr and Ba, in M22. They also 
identified a bimodality among these elements. None of the el¬ 
ements show a correlation with Na, O and Al. The bimodality 
in s-process elements in M22 resembles the case of NGC 1851 
dYong & Grundahj l2008ll although in this case the s-element 
abundance appears to correlate with the Na, Al, and O abun¬ 
dances. The meta l-poor globular cluste r M 15, with [Fe/H] = - 
2.31 + 0.06 dex (ICarretta et al. 2009b), has als o prov en to be 
an interesting case for studv. lSneden et alJ (Il997l[2000h found a 
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scatter of heavy neutron-capture elements but no significant s- 
process enrichment in M15. More recently. iWorlev et al.l (12013h 
found a bimodal distribution of Ba, Eu and possibly La in M15. 
Both modes of the bimodality seem to be indicative of a pollu¬ 
tion scenario dominated by the r-process and only varying due 
to a different degrees of enrichment. A group of La-rich and La- 
poor stars can be identified in our sample but such a small sample 
does not allow us to carryout further investigation. 

The weak i-process accounts for the major fraction of the 
light i-process (Is-) elements, like Y, and occurs in cor e He- 
burning massive stars (M> SMq) (iPignatari et alJ 1201^ . The 
main s-process takes place in thermally pulsating AGB stars 
(lMo<M<8Mo) producing light 5-proc ess (Is-) and heavy s - 
process (hs-) elements, like Ba and La (lArlandini et al.lll999l) . 
Therefore, analysis of ratios between Is- and hs- elements are 
very interesting in order to constrain these formation processes. 
Ligure [14] shows the [/!5/Z5]=[Ba/Y] ratio as a function of the 
heavy 5-process element Ba. Although the large errors should be 
considered, a relation is present between Y and Ba. This rela¬ 
tion could imply a real spread of Y and Ba. Although the spread 
among heavy elements is comparable to the star-to-star scatter 
found in [Le/H], heavy element absorption lines are generally 
difficult to measure. Lor example, Ba absorption lines are quite 
strong among luminous RGB stars so saturation could be a prob¬ 
lem. They have hyperfine structure and involve the choice of at¬ 
mospheric parameters. In particular, they are quite sensitive to 
the choice of microturbulent velocities. These effects could pro¬ 
duce an increase in the uncertainties. 

The 5-process elements can be considered as a signature 
of the processes that occur in intermediate mass AGB stars 
dBusso et al.l2Q0lh . whose wind could have polluted the primor¬ 
dial material from which the second generation of stars formed. 
If the spread in 5-process contents in NGC 4372 is real, this 
would argue against the abundances of 5-process elements be¬ 
ing intrinsic to the cluster and suggest that AGB stars can be a 
possible polluter. However, none of the analyzed heavy elements 
show any clear relations with other lighter elements that could 
suggest an obvious pollution by AGB stars. Since this potential 
spread seems unrelated with the spreads in light elements, the 
variations could be attributed to the original inhomogeneities in 
the gas which formed the cluster. A larger sample survey of stars 
at various evolutionary phases and extensive abundance analysis 
of a range of key nucleosynthetic indicators would be essential 
to confirm whether this potential star-to-star scatter is real. 

5. Summary 

In this paper we present the first detailed chemical abundances 
of 14 elements in 7 red giant members of NGC 4372 using high 
resolution, high S/N spectroscopy. Chemical abundances have 
been computed by Concepcion Node within the GES collabora¬ 
tion. The classical EW method has been used when possible. Lor 
5 elements whose lines are affected by blending, the spectrum- 
synthesis method was preferred. One of the stars of our sample 
shows a radial velocity in agreement with the cluster bulk but 
a distinct chemical signature so we have excluded it from the 
statistical analysis. 

We found a metallicity of [Le/H]= - 2.19 + 0.02 with a cTobs — 
0.03 dex, in good agreement with previous, low-resolution, stud¬ 
ies. We rule out an intrinsic metallicity spread, although the low 
value for the one (excluded) star should be born in mind. We con¬ 
firm the Na-O anti correlation although not very extended, prob¬ 
ably due to our small sample. The abundances of O are relatively 
high compared with other globular clusters which could indicate 


that NGC 4372 was formed in an environment with high O for 
its metallicity. Intrinsic spreads are also seen in other light ele¬ 
ments, in particular an apparent Mg-Al anti correlation has been 
detected. The Le-peak elements generally show good agreement 
with other GCs and halo field stars with no dispersion. The a 
elements show an enhancement of [a;/Le]= -1-0.37 + 0.07 typical 
of other GCs indicating similar fast star formation time scales. 
A relation between light and heavy 5-process elements has been 
identified. A larger sample of stars at various evolutionary phases 
and extensive chemical analysis is required for a more definitive 
analysis. 
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